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ABSTRACT 



Aims. Sloan Digital Sky Survey (SDSS) DR5 photometric data with 120° < a < 270°, 25° < 8 < 70° are searched for new Milky 
Way companions or substructures in the Galactic halo. 

Methods. Five candidates are identified as overdensities of faint stellar sources that have color-magnitude diagrams similar to those of 
known globular clusters or dwarf spheroidal galaxies. The distance to each candidate is estimated by fitting suitable stellar isochrones 
to the color- magnitude diagrams. Geometric properties and absolute magnitudes are roughly measured and used to determine whether 
the candidates are likely dwarf spheroidal galaxies, stellar clusters or tidal debris. 

Results. SDSSJ1000+5730 and SDSSJ1329+2841 are likely faint dwarf galaxy candidates while SDSSJ08 14+5 105, 
SDSSJ0821+5608 and SDSSJ1058+2843 are likely extremely faint globular clusters. 
Conclusions. Follow-up study is needed to confirm these candidates. 
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1. Introduction 



tidal streams in the Milky Way from 
smaller galaxies supports hierarchical 



The discovery of 
the accretion of 
merging cosmogonies (Lynden-Bell & Lynden-Bell 1995] ). 
The large-scale tidal arms of the Sagittarius dwarf 
spheroidal galaxy accretion event were the first to be 
discovered (|Yanny et al. 2003| |Majewski et al. 2003] and 



IRocha-Pinto et al. 2004l i. and more tidal streams continue to be 
detected in Sloan Digital Sky Survey (SDSS; lYork et al. 2000j i 
data. The remarkably strong tidal tails of Palomar 5 are 
discovered by lOdenkirche n et al. 200 II IRockosi et al. 2002I 
and Grillm air & Dionatos 2006a] Another tidal stream 
is connected with N GC5 466 dBelokurov et al. 2006al 
IGrillmair & Johnson 20061 IBelokurov et al. 2006bl re- 
veals the so-called Field of Streams, which shows not 
only Sagittarius arms but also the Virgo overdensity 
dJuric et al. 2005l l. Monoceros Ring (|Newberg et al. 2002|>, 
Orphan stream) Belo kurov et al. 2006bl IGrillmair~2 006a ) 
and a 60°-long stream between Ursa Major and Sextans 
(IGrillmair & Dionatos 20 06b). The fact that tidal streams are 
often associated with globular clusters or dwarf spheroidal 
galaxies implies that it is possible to search for tidal streams by 
identifying new faint satellite dwarf galaxies, globular clusters 
or tidal debris and ascertaining their connections. 

New satellite companions of the Milky Way have been 
discovered in the recent years. About 150 globular clusters 
dHarris 19961 1 and 9 widely-accepted dwarf spheroidal galax- 
ies (dSph) ( Belo kurov et al. 20071 ) were known before the 
SDSS. Since 2005, tremendous progress has been made us- 
ing SDSS data; at least nine new companions were dis- 
covered within two years. Seven of them are dSph satel- 
lites: Ursa Major I ( Will man et al. 2005al l. Canes Venatici 
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I dZucker et al. 2006all, Bootes ([B elokuro v et al. 2006cl Ursa 
Major II (Zuck eret al. 2006bllGril lmair 2006all, ComaBerenics, 
Canes Venatici II, Hercules and Leo IV (Beloku rov et al. 20071 ). 
Two of them are probably globular clusters: Willman 1 
(IWillman et al. 2005bl ) and Segue 1 dBelokurov et al. 20071 ). 

It is impossible to see these objects in SDSS images since 
they are resolved into field stars. They are discovered by de- 
tecting local stellar overdensities in SDSS stellar databases 
(Belokurov et a l. 20071) by adding constraints on magnitude and 
color index. Overdensity study is a quick data mining method to 
find previously unknown satellite galaxies or globular clusters. 
We use a method similar to that of Beloku rov et al. 2006cl but 
utilize different parameters, to scan an area of the sky bounded 
by 120° < a < 270°, 25° < S < 70°, and observe five interest- 
ing overdensities with color-magnitude diagrams similar with a 
globular cluster or a dSph. They are candidate globular clusters 
or dSphs or tidal debris. 

The organization of this paper is as follows. Section 2 de- 
scribes our method of searching for candidates and how we mea- 
sure their properties. Data acquisition and reduction processes, 
including the method for identifying candidates, are discussed 
in section 2. 1 . A template matching algorithm for estimating the 
distance to each overdensity is described in section 2.2. A rough 
measurement of geometric properties and radial profiles are ad- 
dressed in section 2.3. In section 3 we discuss the nature of each 
overdensity. A short conclusion is included in the last session. 

2. Discovery 

2. 1. Data Acquisition and Reduction 

SDSS DR5 provides photometric data in u, g, r, i, and z pass- 
bands (Fukugit a et al. 1996), and c overs 8000 square degrees of 
sky ( jAdelmarTM cCarthy et al. 2007 ). We select only the objects 
with stellar profiles from SDSS databases. SDSS recognizes all 
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point sources as stars, including quasars and faint galaxies. Thus, 
galaxy clusters and arms and halos of bright galaxies contami- 
nate our sample of overdensities. In order to decrease the fre- 
quency of these contaminants, we use only sources with magni- 
tude 19 < i < 22 and color < g— i < 1. We limited the area of 
sky searched to avoid the Sagittarius dwarf spheroidal galaxy's 
tidal arms. Approximately 3.6 x 10 7 point sources were obtained 
in the region 120° < a < 270°, 25° < 5 < 70° from SDSS 
Casjobs[J The magnitudes in each band were corrected for ex- 
tinction using values provided by the SDSS database, which are 



computed following Schlegel et al. 1998 



Data analysis procedures are based on the China Virtual 
Observatory's experimental data access service(VO-DAS). We 
compute the star counts for this data in bins of size 0.2° x 0.2° 
Hence, 750 x 225 bins are generated. Field mean density and 
standard deviation for each bin are defined by the star counts in 
the surrounding bins, within all x 11 window. The standard de- 
viation is derived from all bins in the box except the center bin. 
The strength of the peak in the center bin is estimated from: 



^center — {^center A0/^"j 



(1) 



where n cente r is the number of stars in the center bin, fi is 
the average number of stars per bin in the surrounding bins, and 
a is the standard deviation. For most objects, the bin size is large 
enough to contain an entire star cluster; the surrounding bins 
contain many Galactic stars. For our sample, n cente r is normally 
distributed with fi = 0.33 and a — 0.56. We selected all peaks 
with n cente r > A + 3cr ~ 2. From these statistics, we expect 
about 0.27% of the bins we search will contain statistical fluctu- 
ations high enough to be detected as a peak. We searched 168750 
bins, of which we excluded 18627 zero bins (sky areas which are 
not covered by SDSS). Therefore, we expect about 405 overden- 
sities are actually caused by field star density fluctuation. 

We actually identified 524 overdensity bins, which is 119 
more than would be expected from statistics. The 524 overden- 
sities include statistical fluctuations, previously identified Milky 
Way satellites, bright galaxies, bright stars, and galaxy clusters, 
in addition to the new Milky Way satellites we would like to 
find. We use additional information from the color-magnitude 
distribution of stars in the overdensities to cull out the statis- 
tical fluctuations, as described below. Known objects includ- 
ing bright galaxies, galaxy clusters, and bright stars are man- 
ually identified and eliminated. The known Milky Way satel- 
lites that were re-discovered by this technique include: UMa II, 
UMa I, Willman 1, Pal4, CVn II, CVn I, NGC5272, NGC5466, 
NGC6205, NGC6341, Draco, and Leo III, some of which cover 
more than 2 bins due to their large angular sizes. Overdensity 
bins which are close to the survey boundaries (within 0.5°) are 
also ignored. 

Density contour diagrams and CMDs were generated within 
a 1° x 1° area centered at the remaining overdensity bins. 
Candidates were manually identified by comparing the density 
contour diagram and CMD of each overdensities with those of 
9 known objects: Pal 4, UMa I, UMa II, CVn I, Draco, NGC 
6341, NGC 5466, Willman 1 and NGC 6205. The distance mod- 
ulus (DM) of these templates ranges from 14.2mag to 21.7mag. 
This ensures that the templates cover all possible morphological 
features for a CMD, since the stellar populations that are probed 
within our magnitude limits depend on the distance modulus. 

Five overdensities are finally chosen from over 500 initial 
candidates. Figure[T]displays density contours, core (circle with 



r = 0.1° or r = 0.15°) CMDs, field (annuli area between 
r = 0.5° and r = 0.6°) CMDs and Hess diagrams (subtrac- 
tion between corresponding center area CMDs and field CMDs 
with normalized star counts) for these five candidate Milky Way 
satellites. 



2.2. Distances to the Overdensities 

Distances to the five candidates are estimated by fitting 
isochrone lines to the CMD of each overdensity. An automatic 
template-matching algorithm was developed to discover the 
best-fit Padova isochrone lines(Gi rardi et al. 2004l > for each can- 
didate. Stellar evolution isochrone lines are converted to masks 
where the color index g — i changes from -1 to 2 with step of 
0.03 and the magnitude i changes from -10 to 10 with step of 
0.1. Thus a binary 100 x 200 image mask for each isochrone 
line is formed. The value of each pixel in the image is either 1 (in 
the neighborhood of an isochrone line) or (everywhere else). 
A series of mask templates are created for a range of metallicity, 
age and distance moduli. For each candidate, the template which 
minimizes R, as defined by the following formula, is selected as 
the best match. 



R 



E (C(g-i,i)Q--T(g-i,i,Z,A,DM))) 

g — i,i 

£ C(g-i,i) 

g-%,i 



(2) 
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where C(g — is the candidate's Hess diagram (i = 13 — 
22mag). T(g — i, i, Z, A, DM) is a template mask with specific 
metallicity Z, age A and distance modulus DM. The value of Z 
is 0.0001, 0.0004, 0.001, orO.004. The value of A in log 10 (yr), 
ranges from 9.5 to 10.25 with a step of 0.05. And the value of 
DM ranges from 13mag to 21.9mag with a step of O.lmag. In 
total 5760 templates were evaluated for each overdensity. The 
value of R also measures the goodness of the fit. The smaller R 
is, the better the fit. 

Rmin, the minimum value of R, measures how many stars, 
in percent, are unrelated to the best fit isochrone line in the Hess 
diagrams. In the extreme case that stars are uniformly distributed 
in color-magnitude plane, considering that the isochrone line 
is only a thin line and cannot cover many stars in the plane, 
Rmin ~ 1 is derived. In the opposite extreme case that the 
candidate contains no background stars so that its features are 
completely described by a template, R m in ~ is reasonable. 
Rmin distributions for all known objects, candidates satellites 
and other overdensities are shown in Figure [2] The population 
of known objects and candidates addresses significantly smaller 
R m in values than that of non-candidates overdensities. 

Specifically, we use 9 known globular clusters and dwarf 
spheroidal galaxies: UMa I & II, CVn I & II, Pal4, NGC6205, 
NGC5466, NGC5272, and NGC6341, to test the effectiveness 
of the algorithm. They are all located in the detected area and 
were all detected by our selection process. Their distance mod- 
uli range from 14 to 22, as measured by previous authors (which 
are listed in Table [T]). The standard deviation in the accuracy of 
the DM is a = 0.23mag. 

Best fit isochrone lines are shown together with Hess dia- 
grams of our candidate Milky Way satellites in Figure [3] The 
distance modulus of the best fit template is an estimate of the dis- 
tance to each overdensity. Table|2]displays the parameters of the 
matching templates for each candidates, though we don't expect 
that this method of isochrone fitting will produce accurate mea- 
surements of the ages and metallicities of each candidate. We list 
our best fit metallicity and age parameters simply to show that 



C. Liu et al.: Milky Way Satellites Candidates 



3 



they are not unreasonable for star clusters or dwarf galaxies in 
the Galactic halo. 

Rmin values in table Q] again show the same trend with 
Figure [2] that globular clusters have smaller R m in values, while 
low surface brightness dSphs (for example UMa II) have larger 
Rmin values. The fact that an acceptable DM is derived for UMa 
II, although it has a higher R m i n that reaches the non-candidates 
range(see Figure [2j, gives us confidence that our technique for 
finding distances is valid even for these low surface brightness 
candidates. 



2.3. Properties of Candidates 

Center positions, geometric properties and absolute magnitudes 
are estimated for each candidate and listed in Table [2] 

In order to measure position and geometric properties for 
each candidate Milky Way satellite, stars associated with candi- 
dates are selected. We select stars which are located within the 
1° x 1° around the center position of candidates and also within 
bins in the Hess diagram (right panel in Figure[T]i that have more 
than 0.05 stars per bin. Density contour diagrams of the five can- 
didates are derived by counting selected stars and are displayed 
in Figure [4] 

Because there are relatively few stars in our samples, we did 
not fit for the ellipticity of the candidates. We fit only circularly 
symmetric profiles. Background level is estimated by averaging 
star density of an annulus area around each candidates with ra- 
dius from 30arcmin to 60arcmin. It is then subtracted from radial 
profile before fitting. 

In order to check our methods, we compare the half-light 
radius of Pal5 and Boo using our method and compare the re- 
sults with previous measurements in the literature. For Pal5 
we derive = 2.64' ± 0.07 for the exponential model and 
r/j = 2.91' ± 0.07 for the Plummer model. These values are 
similar to the previously measured value of rh = 2.96' in 
IHarris 19%1 For Boo, we find r h = 14.4' ±1.8 for the expo- 
nential model and = 14.6' ±1.5 for the Plummer model, 
while in Belo kurov et al. 2006cl corresponding values are = 
13.0' ± 0.7 and r h = 12.6' ± 0.7. The results support our 
methods for estimating globular cluster and dwarf galaxy pa- 
rameters, though there may be larger errors for lower luminosity 
objects. All candidates are fitted by exponential and Plummer 
models) McConnachie & Irwin 2006 ). Radial profiles and fitting 
curves are displayed in Figure[5] Half-light radii derived by inte- 
grating exponential and Plummer profiles are tabulated in Table 
[2] Geometrical sizes of these candidates are computed from the 
estimated distances and the angular sizes from each model pro- 
file and also show in Table|2] 

To estimate absolute magnitude of candidates we use Hess 
diagrams derived by subtracting field Hess diagrams from those 
generated by stars inside i'h, normalized by area. Overlapping 
corresponding best matched isochrone masks on these Hess 
diagrams, only star counts located at neighbors of isochrone 
lines are kept. Then, g-band and r-band magnitudes are com- 
puted by integrating all positive flux values in these masked 
Hess diagrams. Total color index (B — V)tot and Total ap- 
parent magnitude Vtot are conducted from Formula [3] and [4] 
( |Fukugita et al. 1996| l. 



(B - V)tot = {gtot - r tot + 0.23)/1.05 
Vtot = rtot + 0.49(B - V) tot - 0.11 
M Vitot = Vtot + 5 - 5log(d) 



(3) 
(4) 
(5) 



Consequently, absolute magnitudes are computed via [5] 
Absolute magnitudes estimated with different radial profile 
models are tabulated in Table [2] As a test, we estimate 
the absolute magnitudes of Com, CVn II, Her and Leo IV, 
which are low luminosity dwarf galaxies discovered in SDSS 
data dBelokurov et al. 20071) . using this method. The maximum 
bias of all four objects between our results and the literature's is 
|AMv,tot| = 0.55mag, which is smaller than the standard er- 
ror(0.6mag) mentioned in the literature. ANOVA analysis also 
shows that the estimation of our method is not significantly dif- 
ferent with those from Beloku rov et al. 20071 



3. Discussion 

It is difficult to identify all 5 candidates reliably using only SDSS 
data; accurate follow up observations using a large telescope is 
required to determine the types, which could be: star clusters, 
dwarf spheroidal galaxies, tidal debris or chance superposition 
of field stars. We expect that these candidates are not merely 
field stars. 

Belo kurov et al. 2007| uses the Rh vs. My plane to identify 
Her, Leo IV, CVn II, Com as dSphs, while Seg I is an unusual 
low luminosity globular cluster. We adopt that same method, 
and plot our candidates with known Milky Way satellites in 
Figure [6] The half-light radii of 5 candidates are derived by 
integrating exponential profiles and Plummer profiles. The 
absolute magnitudes are derived by adding all possible member 
stars flux which located in half-light radius and time 2. We plot 
the 5 candidates in Rh vs. My plane with the mean values 
derived by exponential and Plummer model. Rh and My of 
known globular clusters are come from IHarris 1 996 and those 
of SDSS discovered Milky Way dwarf spheroidal galaxies are 
come from Belokurov et al. 2006c Belok urov et al. 2007 
IZucker et al. 2006bl IWillman et al. 2005al and 

IWillman et al. 2005bl those of dwarf galaxies in the lo- 
cal group are come from llrwan & H atzidimi triou 19951 and 
Mateo 1998 and those of Andromeda dwarf galaxies are come 
from lMcConnachie & Irwin 2006l and lMartin et al. 20061 

SDSSJ0814+5105, SDSSJ082 1+5608 and 

SDSSJ1058+2843: 

SDSSJ08 14+5 105, SDSSJ0821+5608 and 

SDSSJ1058+2843 have surface brightness which are sim- 
ilar to other known satellites discovered from SDSS with 
extremely faint absolute magnitudes, even fainter than AM4, 
the faintest one in Harris 1996 Their location in Figure [6] 
suggests that these stellar systems are more similar to globular 
clusters than to dwarf galaxies. Tidal radii of SDSSJ08 14+5 105 
, SDSSJ082 1+5608 and SDSSJ1058+2843 in Figure [6] are 
computed using the equation 



Ttidal ~ R[ 



ZMmw 



1/3 



(6) 



from |Binney & Tre maine 1987 where R is the candidate's 
galactocentric distance, M can d is its total mass, and Mmw is 
the total mass of the Milky Way within R. We assume that the 
distance from the Sun to the center of the Milky Way is 8kpc and 
v c = 220fcms _1 at distance R. M can d is estimated by compar- 
ing star counts of F and G stars located in main sequences of the 
three candidates and Pal 5 (Figure |7J. Since the formula for the 
tidal radius does not include dark matter, the actual tidal radius 
is larger than estimated if the candidates are dwarf galaxies with 
substantial dark matter content. 
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Notice that SDSSJ08 1 4+5 105 and SDS SJ0821+5608 over- 
lap the Anticenter Stream (Grillmair 2006b), which may be re- 
lated to Monoceros Ring displayed in Figure [8] According to 
|Newberg et al. 2002| the turnoff magnitude of Monoceros Ring 
is at g = 19.4. In Figure [jj we estimate SDSSJ08 14+5 105 
and SDSSJ0821+5608 has a magnitude of the turnoff as g = 
19.5. They are almost at the same distance. However, accord- 
ing to Grillmair 2006b the Anticenter Stream has a turnoff 
magnitude as g = 18.7, quite different from the results of 
Newberg et al. 2002] and substantially closer than the new can- 
didates. The question is, are they globular clusters in the 
Monoceros stream, or are they just part of the debris more com- 
pact than in other areas? 

SDSSJ1000+5730 & SDSSJ1 329+2841: SDSSJ1000+5730 
and SDSSJ1329+2841 are faint, extended objects with half light 
radii similar to that of most of dSphs. They have surface bright- 
ness near \iy ~ Zlmag / 'arcsec 2 , which are lower than all 
known dSphs discovered from SDSS. Their type cannot be de- 
termined only by SDSS due to their faintness, but they look like 
dwarf galaxies. 

4. Conclusions 

In this paper we report five interesting overdensities detected in 
the SDSS database. They show features of globular clusters or 
dSphs in their CMDs. Mass estimation for SDSSJ08 14+5 105, 
SDSSJ0821+5608 andSDSSJ1058+2843 indicate that their tidal 
radii are bigger than their half-light radius, which implies that 
they are likely globular clusters with low surface brightness. The 
half-light radius and absolute magnitude of SDSS J 1329+2841 
suggest a likely dSph. Although surface magnitude measure- 
ments suggest that SDSSJ1000+5730 is fainter than all known 
dSphs, its Hess diagram implies that it is an interesting object 
that needs follow-up observations. 
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Name 


Z from template 


Age from template 

(logi (yr)) 


DM from template 
(mag) 


DM from references 
(mag) 


R mi 71 


References 




UMa II 


0.0004 


10.1 


17.9 


17.5 


0.50 


IZucker et al. 2006b] 


UMal 


0.0001 


10 


19.9 


20 


0.41 


Willman et al. 2005a 


Pal4 


0.001 


10.2 


20.2 


20.02 


0.37 


IHarris 19961 


CVnII 


0.0004 


10.2 


20.9 


20.9 


0.25 


Belokurov et al. 20071 


CVnl 


0.0004 


10.15 


21.7 


21.75 


0.28 


Zucker et al. 2006a 


NGC5272 


0.001 


10.1 


15 


15.04 


0.21 


Harris 1996 




NGC5466 


0.001 


10.05 


16.1 


16.1 


0.34 


Harris 1996 




NGC6205 


0.004 


9.95 


14.7 


14.28 


0.21 


Harris 1996 




NGC6341 


0.001 


10 


14.9 


14.59 


0.28 


IHarris 1996 





Table 2. Parameters of each overdensity. 



Parameter 


SDSSJ0814+5105 


SDSSJ0821+5608 


SDSSJ1000+5730 


SDSSJ1058+2843 


SDSSJ1329+2841 


RA(J2000) 


08 h 13 m 42 s 


08 h 21 m 15 s 


10 h 00 m 28 s 


10 tl 58 m 04 s 


13 tl 29 m 13 3 


DEC(J2000) 


+51°05'27" 


±56°08'16" 


±57° 30' 10" 


±28°42'39" 


+28°41'27" 


l(deg) 


167.743 


161.665 


155.506 


202.649 


45.716 


b(deg) 


33.449 


34.615 


47.372 


64.966 


81.513 


Z 


0.004 


0.004 


0.001 


0.004 


0.0004 


Age(logi (yr)) 


10.5 


10 


10.15 


9.95 


10.1 


(m - M)o(mag) 


15.7 


15.7 


19.6 


16.9 


19.4 


Rmin 


0.5 


0.45 


0.47 


0.36 


0.58 


Distance(kpc) 


ta.st\i 


13.8ll| 


83.2±|| 


9 4 +2.7 


75.9±^ 


rh(exponential)(arcmin) 


6.2 ± 1.0 


4.7 ± 1.0 


8.1 ±2.7 


4.7 ±0.7 
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Fig. 1. Overdensity detections. The first column of diagrams are density contours. Each density contour is constructed by nonpara- 
metric distribution density estimation with a Gaussian kernel and a window width of 0.1°, using stars with 19 < i < 22mag located 
in the 1° x 1° area around the candidate. The second column of diagrams displays the core CMDs of overdensities. Points in this 
column are samples from the smallest circle (r = 0.1° or 0.15°) overlaid on density contours. The third column of diagrams displays 
field CMDs of the overdensities. Points in this column are samples from the annuli between the intermediate circle (r = 0.5°) and 
the biggest circle (r = 0.6°). The right most column of diagrams are Hess diagrams which are the result of subtraction between 
the core CMDs and corresponding field CMDs with normalization on total star counts. The size of each bin in the Hess diagrams is 
0.03mag x O.lmag. 
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Fig. 2. Comparison of R m i n for real objects compared with all original two sigma detections. Two populations of R m i n values 
distribution are displayed. The solid stair line shows all known objects, including globular clusters, dwarf galaxies, bright galaxies 
and galaxy clusters, and our five candidates. The distribution shows two peaks. The peak between 0.2 and 0.3 is due to 5 known 
globular clusters. The peak at 0.4 is associated with dwarf galaxies and normal galaxies. As comparison, the dashed stair line shows 
the distribution of all other overdensities with a gaussian fitting curve(dot line), where the average value is R m in = 0.57 and 
a = 0.1. This figure shows that the globular clusters have the smallest R m in, the dwarf galaxies reach a larger radus, and both 
objects have a significantly smaller R m i n than non-candidates. 
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Fig. 3. The Hess contours of candidates with the best-fitting isochrone lines overlaid on them. Contours are positive densities (0.05, 
0.1, 0.15, 0.2, 0.25, 0.3 stars per bin). The size of each bin is O.OSmag x O.lmag. Best-fitting stellar evolution isochrone lines are 
overlaid with dashed lines. 




Fig. 4. Isodensity contours of five candidates. The contour levels are 2, 3, 4, 5, 6, 7 and 7.5<r above the background. The samples 
that generate the contours are selected from the corresponding Hess diagram. 
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Fig. 5. Density radial profiles for candidates. Background levels are subtracted before fitting models. The solid lines are Plummer 
model fit curves. The dashed lines are exponential model fit curves. 
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Fig. 6. Location of different classes of objects in the plane of half light radius vs. absolute magnitude. Lines of constant sur- 
face brightness are marked. Filled circles with bold labels are the 5 candidates revealed in this paper. Triangles are known 
Globular Clusters dHarris 1996b. Stars with labels are Milky W ay satellites discover ed from SDSS(Belokur ov et al. 20 06c 
IBelokurov et al . 2007 Zucker et al. 2006b, Willman et al. 2005aland lWillman et al. 20"0 5b ). Diamonds are dwarf galaxies in the lo- 
cal groupdrwan & Hatzidimi triou 19951 and IMateo 19 98 ). Rectangles are Andromeda dwarf galaxies(McConnachi e & Irwin 2006 
and lMartin et al. 2006t 
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Fig. 7. Luminosity function comparison between SDSSJ08 14+5 105, SDSSJ0821+5608 and SDSSJ1058+2843, marked as solid 
line, and Pal 5, marked as dot line, within half-light radii. Stars between 4 and 6 mag correspond to spectral types from F to G in 
main sequence. Star counts of F and G stars are used to roughly estimate the total mass. 
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Fig. 8. Overdensities' positions in galactic coordinates in Aitoff projection. The gray area encloses the main area of North Galactic 
Cap covered by the SDSS DR5 data. Filled circles are candidates in this paper. Plus symbols are UMa II, Willman 1, UMa I, CVn 
I, CVn II, Her, Leo IV, Com and Seg 1. Star symbols are globular clusters in the SDSS DR5 area we detect. The triple thin solid 
lines are Anticenter Stream of Grillmair 2006b which may be related to the Monoceros Ring. The thin dotted line is GD-1. The thin 
dashed line is Orphan stream. 



